library(dplyr)

# load cost files
tm_files <- as.list(list.files("data/outputs/", pattern = "k_optimize"))
k_optimize_tm <- lapply(tm_files, function(i) readRDS(paste0("data/outputs/", i))) %>% do.call(rbind, .)

# plot
ldatuning::FindTopicsNumber_plot(k_optimize_tm)
